The role of glutamate in neuronal ion homeostasis: A case study of spreading depolarization.

Simultaneous changes in ion concentrations, glutamate, and cell volume together with exchange of matter between cell network and vasculature are ubiquitous in numerous brain pathologies. A complete understanding of pathological conditions as well as normal brain function, therefore, hinges on elucidating the molecular and cellular pathways involved in these mostly interdependent variations. In this paper, we develop the first computational framework that combines the Hodgkin-Huxley type spiking dynamics, dynamic ion concentrations and glutamate homeostasis, neuronal and astroglial volume changes, and ion exchange with vasculature into a comprehensive model to elucidate the role of glutamate uptake in the dynamics of spreading depolarization (SD)-the electrophysiological event underlying numerous pathologies including migraine, ischemic stroke, aneurysmal subarachnoid hemorrhage, intracerebral hematoma, and trauma. We are particularly interested in investigating the role of glutamate in the duration and termination of SD caused by K+ perfusion and oxygen-glucose deprivation. Our results demonstrate that glutamate signaling plays a key role in the dynamics of SD, and that impaired glutamate uptake leads to recovery failure of neurons from SD. We confirm predictions from our model experimentally by showing that inhibiting astrocytic glutamate uptake using TFB-TBOA nearly quadruples the duration of SD in layers 2-3 of visual cortical slices from juvenile rats. The model equations are either derived purely from first physical principles of electroneutrality, osmosis, and conservation of particles or a combination of these principles and known physiological facts. Accordingly, we claim that our approach can be used as a future guide to investigate the role of glutamate, ion concentrations, and dynamics cell volume in other brain pathologies and normal brain function.

Anions Govern Cell Volume: A Case Study of Relative Astrocytic and Neuronal Swelling in Spreading Depolarization.

Cell volume changes are ubiquitous in normal and pathological activity of the brain. Nevertheless, we know little about the dynamics of cell and tissue swelling, and the differential changes in the volumes of neurons and glia during pathological states such as spreading depolarizations (SD) under ischemic and non-ischemic conditions, and epileptic seizures. By combining the Hodgkin-Huxley type spiking dynamics, dynamic ion concentrations, and simultaneous neuronal and astroglial volume changes into a comprehensive model, we elucidate why glial cells swell more than neurons in SD and the special case of anoxic depolarization (AD), and explore the relative contributions of the two cell types to tissue swelling. Our results demonstrate that anion channels, particularly Cl-, are intrinsically connected to cell swelling and blocking these currents prevents changes in cell volume. The model is based on a simple and physiologically realistic description. We introduce model extensions that are either derived purely from first physical principles of electroneutrality, osmosis, and conservation of particles, or by a phenomenological combination of these principles and known physiological facts. This work provides insights into numerous studies related to neuronal and glial volume changes in SD that otherwise seem contradictory, and is broadly applicable to swelling in other cell types and conditions.

Large extracellular space leads to neuronal susceptibility to ischemic injury in a Na+/K+ pumps-dependent manner.

The extent of anoxic depolarization (AD), the initial electrophysiological event during ischemia, determines the degree of brain region-specific neuronal damage. Neurons in higher brain regions exhibiting nonreversible, strong AD are more susceptible to ischemic injury as compared to cells in lower brain regions that exhibit reversible, weak AD. While the contrasting ADs in different brain regions in response to oxygen-glucose deprivation (OGD) is well established, the mechanism leading to such differences is not clear. Here we use computational modeling to elucidate the mechanism behind the brain region-specific recovery from AD. Our extended Hodgkin-Huxley (HH) framework consisting of neural spiking dynamics, processes of ion accumulation, and ion homeostatic mechanisms unveils that glial-vascular K(+) clearance and Na(+)/K(+)-exchange pumps are key to the cell's recovery from AD. Our phase space analysis reveals that the large extracellular space in the upper brain regions leads to impaired Na(+)/K(+)-exchange pumps so that they function at lower than normal capacity and are unable to bring the cell out of AD after oxygen and glucose is restored.

Dynamics from seconds to hours in Hodgkin-Huxley model with time-dependent ion concentrations and buffer reservoirs.

The classical Hodgkin-Huxley (HH) model neglects the time-dependence of ion concentrations in spiking dynamics. The dynamics is therefore limited to a time scale of milliseconds, which is determined by the membrane capacitance multiplied by the resistance of the ion channels, and by the gating time constants. We study slow dynamics in an extended HH framework that includes time-dependent ion concentrations, pumps, and buffers. Fluxes across the neuronal membrane change intra- and extracellular ion concentrations, whereby the latter can also change through contact to reservoirs in the surroundings. Ion gain and loss of the system is identified as a bifurcation parameter whose essential importance was not realized in earlier studies. Our systematic study of the bifurcation structure and thus the phase space structure helps to understand activation and inhibition of a new excitability in ion homeostasis which emerges in such extended models. Also modulatory mechanisms that regulate the spiking rate can be explained by bifurcations. The dynamics on three distinct slow times scales is determined by the cell volume-to-surface-area ratio and the membrane permeability (seconds), the buffer time constants (tens of seconds), and the slower backward buffering (minutes to hours). The modulatory dynamics and the newly emerging excitable dynamics corresponds to pathological conditions observed in epileptiform burst activity, and spreading depression in migraine aura and stroke, respectively.

Bistable dynamics underlying excitability of ion homeostasis in neuron models.

When neurons fire action potentials, dissipation of free energy is usually not directly considered, because the change in free energy is often negligible compared to the immense reservoir stored in neural transmembrane ion gradients and the long-term energy requirements are met through chemical energy, i.e., metabolism. However, these gradients can temporarily nearly vanish in neurological diseases, such as migraine and stroke, and in traumatic brain injury from concussions to severe injuries. We study biophysical neuron models based on the Hodgkin-Huxley (HH) formalism extended to include time-dependent ion concentrations inside and outside the cell and metabolic energy-driven pumps. We reveal the basic mechanism of a state of free energy-starvation (FES) with bifurcation analyses showing that ion dynamics is for a large range of pump rates bistable without contact to an ion bath. This is interpreted as a threshold reduction of a new fundamental mechanism of ionic excitability that causes a long-lasting but transient FES as observed in pathological states. We can in particular conclude that a coupling of extracellular ion concentrations to a large glial-vascular bath can take a role as an inhibitory mechanism crucial in ion homeostasis, while the Na⁺/K⁺ pumps alone are insufficient to recover from FES. Our results provide the missing link between the HH formalism and activator-inhibitor models that have been successfully used for modeling migraine phenotypes, and therefore will allow us to validate the hypothesis that migraine symptoms are explained by disturbed function in ion channel subunits, Na⁺/K⁺ pumps, and other proteins that regulate ion homeostasis.

Linking a genetic defect in migraine to spreading depression in a computational model.

Familial hemiplegic migraine (FHM) is a rare subtype of migraine with aura. A mutation causing FHM type 3 (FHM3) has been identified in SCN1A encoding the Nav1.1 Na(+) channel. This genetic defect affects the inactivation gate. While the Na(+) tail currents following voltage steps are consistent with both hyperexcitability and hypoexcitability, in this computational study, we investigate functional consequences beyond these isolated events. Our extended Hodgkin-Huxley framework establishes a connection between genotype and cellular phenotype, i.e., the pathophysiological dynamics that spans over multiple time scales and is relevant to migraine with aura. In particular, we investigate the dynamical repertoire from normal spiking (milliseconds) to spreading depression and anoxic depolarization (tens of seconds) and show that FHM3 mutations render gray matter tissue more vulnerable to spreading depression despite opposing effects associated with action potential generation. We conclude that the classification in terms of hypoexcitability vs. hyperexcitability is too simple a scheme. Our mathematical analysis provides further basic insight into also previously discussed criticisms against this scheme based on psychophysical and clinical data.